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We explore the influence of two different polarizations of quantum oxygen vibrations on the spacial 
symmetry of the bound magnetic bipolaron in the context of the t — J model by using exact diag- 
onalization within a limited functional space. Linear as well as quadratic electron phonon coupling 
to transverse polarization stabilize d— wave symmetry. The existence of a magnetic background is 
essential for the formation of a d— wave bipolaron state. With increasing linear electron phonon 
coupling to longitudinal polarization the symmetry of a d-wave bipolaron state changes to a p-wave. 
Bipolaron develops a large anisotropic effective mass. 

PACS numbers: 71.27.+a,71.38.Mx, 71.38.-k,74.20.Rp 



Soon after the discovery of high-T c superconductivity 
the quest for the pairing mechanism focused on magnetic 
fluctuations due to a broadly accepted conjecture that 
phonon mechanism alone is not strong enough to pro- 
duce high transition temperatures as observed in high- 
Tc compounds. Recently, a growing evidence is emerg- 
ing in favor of the significance of lattice degrees of free- 
dom in high-T c compounds [H, 0, 01 . The interplay be- 
tween strong correlations and lattice degrees of freedom 
[3| seems to be responsible for many unusual properties 
of cuprates in the low doping regime, such as kinks [H, 
stripes @, and waterfall [EI3 structures. 

A long-standning objection against phonon-based 
mechanism for high-T c superconductivity is based on a 
widely accepted notion that coupling to phonon degrees 
of freedom is predominantly consistent with s— wave pair- 
ing, not characteristic for cuprates. Despite recent dis- 
covery that weak EP coupling to acoustic phonons in the 
presence of large on-site Coulomb interaction leads to 
d— wave pairing [8|, the role of short-wavelength oxygen 
oscillations on the symmetry of the paired state remains 
a challenging open problem. 

Investigation of correlated models coupled to phonons 
were based on exact diagonalization (ED) calculations 
on small lattice systems [3, E3, EH, E3, slave-boson ap- 
roache s 1131 1 14 . Il5 1 , dynamical mean-field calculations 
LZI, UM, coherent states Lanczos method UM, and 
quantum Monte Carlo methods (QMC) [H, El. In 
Ref. authors present a detailed study of the influ- 
ence of EP coupling to a single phonon mode on for- 
mation of inhomogeneous charge structures in the t — J 
model. They show that half-breathing mode stabilizes 
a stripe phase. In contrast, using slave-boson approach 
authors of Ref. 13!] suggest, that half-breathing mode en- 
hances d— wave pairing. They furthermore underline 
the importance of off-diagonal EP coupling modulating 



the hopping and the spin-exchange terms. In contrast, 
authors of Ref. [llj find that diagonal EP terms exceed 
off-diagonal ones by nearly two orders of magnitude. 

ED calculations of the t — J model show that the 
d— wave symmetry of a bipolaron [22 
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is not robust against addition of longer range hopping 
terms 27] while recent QMC calculations of the Hub- 
bard model yield T c far below those of cuprates [28j |. 
There seems to be a need to uncover additional mech- 
anism that would help stabilize the d— wave symmetry of 
a bound bipolaron state. In this Letter we show that EP 
coupling to a transverse polarization (TP) of oxygen vi- 
bration provides an important mechanism that stabilizes 
the d— wave symmetry of a bound hole pair with a small 
effective mass. 

We solve a system of two holes in the t — J model de- 
fined on an infinite two-dimensional lattice by extending 
the method for a single hole based on exact diagonaliza- 
tion within a limited functional space 0, [13| . We intro- 
duce diagonal EP coupling to either TP of oxygen (O) 
vibration relevant for the description of buckling modes 
or longitudinal polarization (LP) of O vibration relevant 
for description of bond-streching modes. We investigate 
the following Hamiltonian 
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,) is a projected fermion opera- 



tor, t represents nearest neighbor overlap integral, the 
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sum runs over pairs of nearest neighbors, ai are 

phonon annihilation operators and ni = X)s n M- The 
third term represents linear EP coupling to LP of O vi- 
bration with respect to Cu-O-Cu bond, see also Fig.QJa). 
Fourth term is chosen cither linear (/3 = 1) or quadratic 
(}3 = 2) in O displacement, describing TP of O vibration, 
Fig. [TJa) . j3 — 2 is chosen to describe the CuO plane with 
no pre-buckling of O positions. Sums over 5 in the latter 
two terms run over two orthogonal nearest neighbor Cu 
positions. Lattice vibrations on O sites are independent 
- we do not predispose any particular phonon mode with 
the exception of limiting our calculation to either TP or 
LP of O oscillation. In treating quantum phonons we 
follow well established approach of Ref. [30( ■ 

The construction of the functional space starts from a 
Neel state with two holes located on neighboring Cu sites 
and with zero phonon quanta. Such a state represents 
a parent state of a translationally invariant state with 
a given momentum k. In the case of a high symmetry 
point k = (0, 0) the parent state can be chosen to have 
d— , s—, or p— wave symmetry as for the case of d— 
and s— shown in Fig.[TJb). The starting state is written 
as \4>^) a = J2-y(~ l) Ma '" ) '^coc- ) .|Neel; 0), where sum runs 
over four nearest neighbors in the case of d— and s— wave 
symmetry and over two in the case of p x ( y yw&ve while 
Mail), a G {d,s,p} sets the appropriate sign. 

We generate new parent states by apply- 
ing the generator of states |l0;' l '*' ) )a} = 
H kin + Hj + H ph y h |0(°)) a ; n h = l,...,N h where 

#kin represent the first term in Eq. [1] Hj denotes a 
part of the second term in Eq. [T] which is only applied to 
erase spin flips that were generated through succeeding 
application of ifkin, as for a particular case depicted in 
Fig. [He). iJph represents either third or fourth term in 
Eq. [1] This procedure generates exponentially growing 
basis of states, consisting of different shapes of strings in 
the vicinity of the hole with maximum lengths given by 
as well as phonon quanta created along paths of both 
holes. Identical basis functions, generated by different 
processes, are chosen only once. We have used Nh = 8 
that lead to N st = 13 x 10 6 states. Full Hamiltonian in 
Eq. [1] is diagonalized within this limited functional space 
taking explicitly into account translational symmetry. 

In Fig. [2] we present the energy difference E p — Ed 
between the lowest p— and the d— wave state for two 
different values of J/t as a function of qp/t and g/t for 
the case of TP and LP respectively. At J/t = 0.1 the 
two-holes are unbound at qi^ = g = [H, [23| and de- 
generate p— wave ground state is found, E p — Ed < 0, 
see Figs. [2ja) and (b). Increasing qp/t and g/t leads 
to rather surprisingly distinct results. In both cases in- 
creasing EP coupling leads to a formation of a bipolaron, 
as also evident from Figs. Wifi) and (f) and the discus- 
sion later in the text. While coupling to TP leads to a 
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Figure 1: (Color online) (a) Schematic representation of LP 
and TP vibrations of O atom (middle) with respect to Cu- 
O-Cu bond, (b) schematic representation of a d x 2_ y 2 — wave 
(s— wave) top (bottom) signs two-hole starting wavefunction. 
Fermion sign convention places the first hole depicted with the 
full circle to the left-most, if the pair is vertical, then bottom- 
most position. Only Cu sites are presented with small dots in 
(b); (c) schematic representation of succeeding applications 
of different off-diagonal parts of Hamiltonian in Eq. [1] for 
rih = 4, starting from a single hole-pair in the Neel state with 
zero phonon quanta. Dots represent Cu and O atoms, holes 
are denoted by open circles, crosses represent spin-flips (over- 
turned spins with respect to the original Neel configuration 
of spins, localized on Cu sites), vertical arrows indicate and 
point to the numbers of excited TP phonon quanta. 



formation of a bound state with the d— wave symmetry, 
coupling to LP in contrast favors a bound state with the 
p— wave symmetry. This effect is even more pronounced 
at larger value of J/t = 0.4 where at 91,2 = Q = a bound 
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magnetic bipolaron is already formed [2j 
d— wave symmetry. By increasing qi/t, d— wave symme- 
try is stabilized. Increasing linear EP coupling q\ /t leads 
to an initial increase of E p — Ed followed by a decrease, 
E p — Ed — ► around q\/t > 0.75 due to a crossover to a 
strong EP coupling regime, Fig. [2](a) . In contrast, linear 
EP coupling to LP drives even a bound rf— wave bipo- 
laron state at J/t = 0.4 and g = to a bound state with 
a p— wave symmetry at g/t ~ 0.61, see Fig. HJb). 

Effective bipolaron mass m aa = t [9 2 -E(k)/9k9k] ^, 
computed in its eigen-directions, presented in Figs. |^c) 
and (d), is isotropic in the case of d— wave symmetry and 
anisotropic, with the anisotropy ratio m yy /m xx ^3 — 10 
in the case of p— wave state. At J/t = 0.4 m aa further- 
more shows only a weak increase with q-z/t, see Fig.^c). 
Even more surprising is the decrease of the effective mass 
at J/t = 0.1 in the regime of a bound wave bipolaron, 
i.e. for q^/t > 0.5. Note, that the nonanalytic behavior 
of m xx is a consequence of the symmetry change from p— 
to d~ state at q2/t ~ 0.5 as also seen from Fig. (2fa). 

Focusing on linear EP coupling to LP, m xx at J/t = 0.4 
starts a rapid increase signaling the approach to strong 
EP coupling regime just below the transition to the 
p— wave state, around g/t ~ 0.57, see Fig. [^d). As 
the system enters p— wave state the mass again becomes 
anisotropic. 

In Fig. [3] we present the probability of find- 
ing a hole-pair at a distance of r: P{r) = 
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Figure 2: (Color online) E p — Ed at tJo/i = 0.2, and two 
different strengths of J/t vs. qp/t in the case of EP coupling 
to TP (a) and vs. g/t in the case of linear EP coupling to LP 
(b). Effective masses m aa vs. q% /t in (c) and vs. g/t in (d). 
The ground state wave-vector is k = except in the regime 
g/t> 0.57 where k = (tt, 0). 



(E(i^j) n i l «j i<5 [l i -j|- r ])/(E(i^j)«Nj l > ) and average 
hole distance (d) = ^ r rP(r). We first focus on the effect 
of EP coupling to TP, see Figs. [SK^e). At J/t = 0.1, 
bipolaron is unbound in the regime (q 2 ft < 0.5, and 
qi/t < 0.22), nevertheless, (d) remains finite due to a 
limited Hilbert space where the maximal inter-hole dis- 
tance is given by Z max = Nh + 1 — 9. Increasing Nh would 
lead to further increase of (d) in this regime as well as 
to further spread of P(r) towards larger r, see Fig. [3fc) 
for J/t = 0.1 and qi j2 = 0. In this range of parameters 
we observe no exponential decay of P(r), see Fig. [3je). 
In contrast, in the regime of a bound bipolaron, i.e. for 
J/t = 0.1 and (q 2 /t > 0.5 and q 1 /t > 0.22) as well as at 
J/t = 0.4, (d) and P(r) do not change much with further 
increasing Nh and exponential decay is clearly observed 
in Fig. [IJe). Our tests performed on smaller systems 
(Nh = 4 and 6) reaffirm that results in the regime of a 
bound bipolaron have indeed converged close to a ther- 
modynamic limit. Good agreement of P(r) at J/t = 0.4 
and qig = g = is found with ED calculation on 32-sites 
cluster, Ref. [22|,[24|]. Structure of a bound bipolaron, re- 
vealed by P(r) at J/t — 0.4 and qi 2 = 0, is remarkably 
similar to that computed at J/t = 0.1 and q 2 /t = 1.0 
and (3 = 2, Fig. EJc). Both, quadratic and linear EP 
coupling to TP lead to a formation of a bound bipolaron 
with the d— wave symmetry even in the case of small 
exchange interaction J/t = 0.1 where magnetic mech- 
anism is not strong enough to form a bound magnetic 
bipolaron. To investigate whether coupling to TP alone 
can lead to d— wave state in the absence of a magnetic 
background, we have solved a problem with two spinlcss 
particles quadratically coupled to TP or linearly to LP 
using topology of a Cu-0 plane. By increasing q-i/t or 
g/t we obtain in both cases a bipolaron with a p— wave 
symmetry, Figs.[3ja) and (b). We thus emphasize an im- 
portant conclusion: EP coupling to TP stabilizes d— wave 
symmetry of a hole-pair, however, the existence of a mag- 



netic background as found in the t — J model seems to 
be essential precondition for the formation of a d— wave 
state. 
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Figure 3: (Color online) Average hole distance (d) computed 
at luo /t — 0.2 and two different values of J/t vs. qp/t;/! — 1,2 
in (a) and vs. g/t in (b). Results for the system of spinless 
particles (SL) are shown with dotted lines in (a) and (b); P(r) 
at chosen values of q\,i/t in (c) and g/t in (d). P(r) is nor- 
malized to 5Z r P(r) — 1; corresponding exponential scalings 
of P(r) for r > 2 are shown in (e) and (f). 

Turning to linear EP coupling to LP phonons we start 
from J/t = 0.1. A p— state of two separate holes changes 
to a bound bipolaron state with increasing of g/t at about 
g/t ~ 0.5, see Fig. [3jb). Transition to a bound state is 
not sharp as in Fig.[3Ja) since there is no change of a sym- 
metry. Nevertheless, in the Nh — > oo limit, we anticipate 
a sharp transition from an unbound to a bound bipolaron 
state. Stabilization of a p— wave state under the influence 
of LP is even more evident when starting from a d— wave 
bound bipolaron state at J/t — 0.4. With increasing 
g/t, a change of symmetry occurs around g/t ~ 0.61, 
from a d— to a p— wave state, see Figs. [3Ub,d,f). A de- 
tailed inspection of a bound p— wave state in the regime 
g/t > 0.61 reveals unusually simple structure where the 
probability of finding holes at a distance r = 2 is more 
than 0.6. This is in a sharp contrast with the structure 
of a d— wave bound state where P(r = 2) < 0.1 and 
the maximal value of P(r) is at r — compare also 
Figs. BJa) and (b). Our calculations of hopping term 
modulated by LP phonons, as suggested in Ref. [13j, as 
well leads to stabilization of a bipolaron with a p-wave 
symmetry, nonetheless, with a distinct spacial structure. 

From Fig. [3] is as well evident that linear coupling 
to TP leads a stronger attraction between holes than 
coupling to LP as seen from Figs. [3te) and (f) that 
show steeper decay of P(r) for TP at comparable values 
qi/t = g/t = 0.7. At small J/t = 0.1 a bound bipolaron 
state is obtained at unexpectedly small value of the di- 
mensionless EP coupling constant X q = qj/4ujot ~ 0.06 
in the TP case in contrast to A g = g 2 /Auoot ~ 0.31 in the 
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LP case. This result suggests that even a small pre- 
buckling within the CuO plane, that generates non-zero 
linear EP coupling term to TP, may have a pronounced 
effect on the attraction between magnetic polarons. 
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Figure 4: (Color online) (a) C(r) and Af(r) at J/t = 0.4, 
qijt = 1.0, and /3 = 2 (d— wave). Radii of circles, representing 
C(r), located on Cu-sites are, proportional to the probability 
of finding a hole-pair at a distance of r. Empty circle is located 
at r = 0. Sides of squares, representing 7V(r), located on O- 
sites, are proportional to average numbers of phonon quanta 
at a distance r from a hole; (b) C(r) and X(r) at J/t = 0.4 and 
g/t = 0.7 (p— wave). Lengths of arrows, representing X(r), 
are proportional to displacements of O atoms along Cu-Cu 
bonds relative to the hole position. 

To investigate in more detail the nature of 
the magnetic-lattice bipolaron, we simultaneously 
present two correlation functions: hole- hole density 
C(r) = (X^i n i n i+r)i an d hole-phonon number Af(r) = 

(Ei^Ni+r^+r)' m Fig- Ufa), f° r t ne case °f a bound 
d— wave bipolaron state. Largest phonon numbers are 
found at the closest possible distance from the hole. The 
structure of C(r) is consistent with d x 2_ y 2 symmetry de- 
spite its largest value at a distance of r = as already 

In Fig. IUh) we show C(r) 

<2j +r ) ) , measuring displace- 
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pointed out in Refs 

and X(r) = (Ei^(4+ 
ments along Cu-Cu bonds relative to the position of the 
hole, for the case of a p— wave ground state. Both cor- 
relations display the unidirectional spacial distribution. 
Correlation functions, presented in Figs. Ufa) and (b), 
show detectable values only up to r < 3, despite maximal 
distance l max = Nh + 1 = 9, allowed in our calculations. 

In conclusion significantly different bipolaron states 
are found when EP coupling to either TP or LP is 
switched on. Linear as well as quadratic EP coupling to 
TP stabilizes a d— wave bipolaron state. The magnetic 
background is essential for the formation of a d— wave 
bipolaron. The effective bipolaron mass remains small in 
the case of quadratic EP coupling despite lattice driven 
binding of the bipolaron. 

In contrast, increasing linear EP coupling to LP 
phonons changes the symmetry of a bound bipolaron 
from a d— wave state at zero EP coupling to a p— wave 
state followed by a substantial change of the density- 
density correlation function. Since this state also has 
a large and anisotropic effective mass and unidirectional 
spacial distribution we may speculate, that in a system 



with finite doping linear EP coupling to LP of O vibra- 
tion would lead to formation of stripe states. This finding 
is consistent with inelastic neutron experiments show- 
ing strong coupling to the bond-streching mode in and 
around the vicinity of the stripe phase in copper oxide 
superconductors |5j. 
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